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EXECUTIVE  SUMMARY 


OBJECTIVE 

Estimate  the  array  gain  achievable  using  adaptive  matched-field  beamforming  of  large  multi¬ 
dimensional  acoustic  arrays  in  coherent  and  incoherent  noise  Eelds. 

APPROACH 

A  mathematical  description  of  adaptive  matched-field  beamforming  for  coherent  and  inco¬ 
herent  noise  fields  is  developed.  Data  covariance  matrices  are  defined  for  150  interference  noise 
sources  varying  in  level.  The  matched-field  beamformer  output  and  array  gain  are  estimated 
when  the  array  is  focused  on  each  of  the  sources  for  various  interference-to-independent-noise 
ratios. 

RESULTS 

Array  gain  is  predicted  to  be  in  excess  of  directivity  index  by  the  interference-to- 
independent-noise  ratio  and  is  observable  for  ratios  of  0  dB  and  greater. 

RECOMMENDATIONS 

Experiments  using  multidimensional  arrays  of  100  to  1000  hydrophones  are  recommended  to 
test  the  predictions. 
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INTRODUCTION 

Matched-field  beamforming  (MFB)  is  a  generalization  of  plane-wave  beamforming  in  which 
the  complex  steering  vectors  are  modeled  as  solutions  to  the  wave  equation  for  the  measured 
environment  rather  than  as  plane  waves.  Both  linear,  or  Bartlett,  and  adaptive  (e.g.,  minimum 
variance  (MV))  processing  schemes  have  been  reported  in  the  literature.^’^^>^*^>^  Beamformers 
provide  signal-to-noise  improvement  (i.e.,  array  gain)  and  estimates  of  signal  parameters.  In  con¬ 
trast  to  plane  wave  beamforming  of  a  volumetric  array,  which  estimates  the  azimuth  and  eleva¬ 
tion  of  the  signal  arrival,  MFB  estimates  the  range,  depth,  and  azimuth  of  the  source  of  signals 
incident  on  the  array. 

This  report  reviews  linear  and  adaptive  matched-field  beamforming  processors  and  derives 
expressions  for  the  mean  array  gain  of  the  processors,  assuming  either  incoherent  or  correlated 
noise.  The  energy  contribution  of  an  interference  on  a  sidelobe  is  calculated  to  provide  peak 
sidelobe  criteria  for  anay  design.  The  sensitivity  to  mismatch  of  the  processors  is  examined  ver¬ 
sus  the  number  of  hydrophones  and  the  input-signal-to-noise  ratio. 


MATCHED-FIELD  BEAMFORMERS 

Consider  a  multidimensional  array  with  horizontal  and  vertical  aperture  sufficient  to  resolve 
interferences  and  element  spacing  larger  than  V2  in  an  arbitrary  ocean  environment  similar  to 
that  shown  in  figure  1.  The  acoustic  array  data  contain  samples  of  the  signal  and  noise  incident 
upon  the  array.  The  hydrophone  data  are  conditioned,  filtered,  sampled,  and  Fourier  analyzed, 
generating  a  complex  Fourier  coefficient  for  each  frequency,  hydrophone,  and  time  sample.  The 
data  for  each  frequency  bin  are  organized  into  a  data  matrix; 
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where  N  is  the  number  of  array  hydrophones  and  T  is  the  number  of  time  samples.  A  covariance 
matrix  is  estimated  by  averaging  over  the  time  that  targets  are  in  a  search  cell  (■•■  denotes  com¬ 
plex  transpose): 

K  =  XX+  (2) 


The  signal  energy  arriving  during  the  averaging  time  in  a  search  cell  is  estimated  using  either 
linear  or  adaptive  matched-field  beamforming. 
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Figure  1.  Matched-field  beamforming  process. 


UNEAR  MATCHED-FIELD  BEAMFORMING 

Linear  beamforming  in  the  frequency  domain  estimates  the  received  energy  by  applying 
steering  vectors  to  the  estimated  covariance: 

{  (a))  =  s-^(a)  k  s(a)  (3) 

where  s'*"  (a )  is  the  complex  steering  vector.  The  energy  is  estimated  for  each  search  cell  in  the 
search  area  by  varying  the  trial  parameter,  a ,  of  the  steering  vector  and  searching  for  an  energy 
peak.  The  position  of  the  peak  is  an  estimate  of  the  range,  depth,  and  azimuth  of  the  signal 
source.  s(a  )  is  calculated  by  using  measurements  of  the  array  element  locations,  estimates  of  the 
sound  speed  structure  derived  from  oceanographic  measurements,  and  acoustic  propagation 
models. 

In  the  far-field  plane- wave  beamforming  case,  the  signal  model  assumes  a  plane  wave  para¬ 
meterized  in  azimuth  and  elevation,  a  —  (0, 4>).  For  focused  beamforming,  the  signal  model  is  a 
spherical  wave  propagating  from  a  target  location  parameterized  in  range,  azimuth,  and  eleva¬ 
tion,  a  =  (r,  0,  (j)).  In  both  cases,  the  complex  steering  weight  is  of  the  form  Sj  =  e'i^j  with  unit 
amplitude  and  is  normalized  according  to 


=  N  (4) 

;=i 

For  matched-field  beamforming,  the  signal  is  modeled  with  an  acoustic  propagation  model 
parameterized  in  trial  range,  depth,  and  azimuth,  c  =  r  =  (r,d,  0)  is  the  set  of  search  space 
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parameters.  The  vector  s  =  ( si,  S2>  •  •  ■W )  is  called  the  replica  for  a  signal  at  the  trial  location. 
Using  the  normalization  from  Eq.  4,  the  steering  vector  is  derived  from  the  modelled  pressures 
by 


pM) 


(5) 


The  replica  accounts  for  the  propagation  multipaths  and  curved  wavefronts  included  in  the 
propagation  model.  The  replica  amplitude  and  phase  both  depend  on  the  hydrophone,;,  because 
of  Ae  depth  dependence  of  the  propagation  and  vertical  multipath  interference  pattern.  Plane- 
wave  and  focused  beamforming  are  special  cases  of  matched-field  beamforming  which  assume 
plane-  and  spherical-wave  propagation  models. 

ADAPTIVE  MATCHED-FIELD  BEAMFORMING 


In  adaptive  MFB,  the  output  energy  is  estimated  according  to 

(  =  w  +  (n)  k  w(fl) 


(6) 


A 

The  adaptive  weights  are  determined  from  the  estimated  covariance,  K,  to  minimize  the  out¬ 
put  energy  subject  to  the  constraint 

N 

w+(a)  s(fl)  =  Sjia)  =  1  (7) 

Adaptive  MFB  minimizes  the  energy  received  by  nulling  interferences,  while  the  constraint 
assures  that  signals  matching  the  replica  vector,  s,  are  passed  without  distortion.  The  minimum 
variance  weight  vector  is  given  by 


= 


k  ^(a) 

s  +  (n)k~^(o) 


and  the  output  energy  is  obtained  by  substituting  in  Eq.  6: 

1 


(  >  = 


-1 


s  +  (a)  K  s(a) 


(8) 


(9) 


Minimizing  energy  is  accomplished  in  the  process  by  using  the  measured  data  to  estimate  the 
amplitudes  and  phases  of  interference  signals  which  do  not  match  the  steering  vector  and 
coherently  subtracting  the  energy.  Knowledge  of  the  replica  vectors  for  the  interferences,  includ¬ 
ing  the  associated  locations  and  propagation  parameters,  is  not  required  by  the  process.  Knowl¬ 
edge  of  the  interference  positions  can  be  used  in  the  inverse  problem  to  estimate  the  sound-speed 
structure  and  array  element  positions.^ 
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SUMMARY 

Matched-field  beamforming  allows  the  use  of  large  vertical  apertures  which  fill  the  water 
column.  It  provides  the  potential  of  separating  submerged  target  signals  from  surface  and  sound 
channel  noise.  If  plane-wave  beamforming  were  to  be  performed  on  the  array,  its  vertical  resolu¬ 
tion  would  be  ~1  degree.  Arrival  angles  from  nonbottom  interacting  sources  typically  vary 
between  ±20  degrees.  The  plane-wave  beamformer  smears  the  signal  multipath  arrivals,  result¬ 
ing  in  gain  degradation.  MFB  combines  the  coherent  vertical  signal  multipath  and  corrects  for 
changes  of  signal  arrival  angle  over  the  vertical  aperture. 

Adaptive  MFB  reduces  sidelobe  energy  by  minimizing  the  energy  that  does  not  match  the 
signal  replica  constraint.  By  focusing  on  the  signal  source  and  minimizing  the  sidelobe  energy, 
the  processor  increases  the  signal-to-noise  ratio  (SNR). 


ARRAY  GAIN  IN  NONUNIFORM  INCOHERENT  NOISE 

For  the  purposes  of  this  section,  it  is  assumed  that  the  ocean  has  a  single  point  source  with 
incoherent  noise  which  varies  from  element  to  element.  K  then  has  the  form 


K  =  Ojdd'*’  +  Off  A 


(10) 


where  dd"*"  is  the  dyadic  matrix  formed  from  the  vector  of  pressures  from  a  source  at  the  source 
position  r,  =  (r„  d„  6^.  The  average  signal  energy  incident  on  the  array  is  O5 .  A  is  a  diagonal 
matrix  of  the  relative  noise  energy  and  o^  is  the  average  noise  in  the  array  elements.  If  the  noise 
energy  is  equal  on  all  of  the  hydrophones,  A  =  I,  the  identity  matrix. 

LINEAR  MFB  IN  INCOHERENT  NOISE 


For  a  linear  beamformer,  the  output  energy  from  Eq.  3  is  given  by 

(Ei(a,  f,))  =  of  s'^(a)d(r5)d  +  (r,)s(a)  -I-  of  s'^ (a)  As(fl) 
Assuming  perfect  match,  s  =  d,  and  using  Eq.  4,  the  output  energy  becomes 

{Ei(,a,7s))  =  N^of  -1-  N 

where 

N 


r-s  _  S  +  As  ;  =  1 


N 


N 


(11) 


(12) 


(13) 


is  the  signal  weighted  average  relative  noise  power  in  the  array  elements,  which  accounts  for  the 
combination  of  noise  and  signal  variation  across  the  array.  For  example,  if  the  noise  is  high  on 

— S 

the  elements  where  the  signal  energy  is  high,  then  Xn  will  be  higher  and  the  gain  will  be  lower. 
Note  that  if  either  the  noise  is  equal  on  all  elements  or  all  the  signal  powers  are  equal  (as  in 

—S 

plane-wave  beamforming),  then  X„  =1. 
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From  Eq.  12,  we  see  that  when  the  matched-field  beamformer  is  steered  at  the  signal,  the  sig¬ 
nal  energy  increases  by  and  the  noise  energy  increases  by  N.  The  array  gain  is  defined  as  the 
output  SNR  divided  tjy  the  input  SNR; 


AG  =  10  log  - 


/<£(oi  =  0))\ 
\<£(o?  =  0)); 


For  linear  MFB,  the  array  gain  is  given  by 


AGi  =  10  log 


te) 


The  array  gain  in  uniform  incoherent  noise,  A  =  I,  is  the  directivity  index,  £>/  =  10  log  N. 
ADAPTIVE  MFB  IN  INCOHERENT  NOISE 

For  adaptive  MFB,  the  output  energy  from  Eq.  9  is  given  by 

{£^(a,r,))  = - ! - ;; - 

s+(d)  (o|  d(r.)d  +  (n)  +  oJA)  s(i) 

Assuming  perfect  match  (s  =  d),  and  using^ 

[usv+  +  b]'^  =  B-‘-B-‘u[v+B-iU  -I-  S-1]"*V+B-* 

the  output  energy  and  optimum  weight  vector  become 


{E^ia,rs))  =  of  +  ^ 


A~^s 
s  +  A'^s 


r-TT'*  S‘'‘A“^S  j  =  l 

where  K„  = - —  is  the  mean  of  the  inverses  of  the  relative  noise  power  in  the 

array  elements  weighted  by  the  relative  signal  powers.  The  array  gain  becomes 

AG^  =  10  log  (20) 


If  the  noise  is  equal  on  all  elements,  X„  '  =1,  adaptive  MFB  has  the  same  weights  and  array 
gain  as  linear  MFB.  When  the  noise  varies  from  element  to  element,  the  elements  with  high 
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noise  are  weighted  less  than  the  elements  with  tow  noise.  Thus,  even  in  the  case  with  incoherent 
but  unequal  noise,  the  array  gain  is  increased  with  adaptive  processing. 

ARRAY  GAIN  IN  NOISE  WITH  MULTIPLE  UNCORRELATED 

INTERFERENCES 

In  this  section,  the  array  gain  is  derived  by  assuming  that  the  array  noise  consists  of  signals 
from  M  uncorrelated  sources  plus  uniform  independent  element  noise.  The  source,  labeled  k,  is 
defined  as  the  signal  upon  which  the  array  is  focused.  The  covariance  matrix  for  this  assumption 
is 

K  =  +  K„ 

M 

K,  =  +  oil  (2i; 

>=1 
j^k 

It  is  also  assumed  that  the  independent  noise  is  equal  on  all  elements.  The  effects  of  unequal  but 
independent  noise  are  similar  to  that  discussed  in  the  preceding  section. 

LINEAR  MFB  WITH  MULTIPLE  INTERFERENCES 

Assuming  that  the  array  is  steered  at  the  signal  s  =  djt,  the  energy  at  the  output  of  a  linear 
matched-field  beamformer  (using  Eq.  3, 4,  and  21)  is  given  by 

{£i)  =  JV2oJ  +  JV2|;oJi^,  +  AroS  pj, 

j^k 

where  s'djd,s  =  N^Lj,  andLyj  is  the  sidelobe  power  level  for  a  linear  beamformer  of  the  signal 
when  the  array  is  steered  at  s. 

For  this  case,  the  array  gain  is  defined  as  the  output-signal-to-total-noise  ratio  divided  by  the 
input-signal-to-total  noise-ratio 


AG  =  10  log 


/(£(a2  =  o;  =  0))\ 

[  (£(o^  =  0)>  j 


where  al  =  ol+  ^  oj  is  the  total  noise  energy.  The  linear  array  gain  is  given  by 


')  +  N 


AGj^  =10  log  TV-IO  log 


1  +  aj 


6 


TR1661 


where  af=  af/ai  is  the  input  interference-to-independent-noise  ratio  (IINR)  for  the  interfer 

r  ^ 


ence  and  10  logu/  =  lO  log 


I-/* 


is  the  total  IINK,  IINRt. 


Note  that  array  gain  is  reduced  from  the  isotropic  noise  case  by  the  interference  signals  in  the 
sidelobes.  If  there  is  no  interference  in  the  main  beam  and  all  the  sidelobes  are  lower  than  1/N, 
gain  better  than  10  log  N  can  be  achieved  for  all  resolved  signals.  When  the  anay  is  steered 
where  an  interference  is  in  a  high  side  lobe,  the  gain  is  less  than  10  log  7^.  Obtaining  low  side¬ 
lobes  requires  a  large  number  of  elements,  closely  spaced  and  with  high  tolerance  on  the  ampli¬ 
tude  and  phase.  Linear  MFB  gain,  even  with  most  careful  array  design,  will  typically  be  less  than 
10  log  N  because  sidelobes  within  a  few  decibels  of  the  main  lobe  occur  at  convergence  zone 
(CZ)  spacing.^  These  sidelobes  result  from  propagation  in  the  oceanic  waveguide  and  cannot  be 
reduced  by  array  design  methods. 

In  the  discussion  of  linear  mismatch  effects,  Eq.  42,  it  is  shown  that  the  mean  sidelobe  level 
is  17  =  ^  .  The  linear  MFB  array  gain  in  coherent  noise,  assuming  all  interferences  are  at  the 
mean  sidelobe  level,  reduces  to  the  directivity  index,  aUI  =  lO  logjv  =  Di. 

ADAPTIVE  MFB  WITH  MULTIPLE  INTERFERENCES 

Adaptive  MFB  reduces  the  sidelobe  energy  whenever  the  sidelobe  is  greater  than  a  few  deci¬ 
bels  down  from  the  main  lobe.  This  helps  in  removing  the  CZ  sidelobes  and  allows  the  use  of 
thinned  arrays  with  element  spacing  greater  than  X/2.® 

Using  Eq.  17,  it  can  be  shown  that  if  K  =  a  dd+  +  Kn.  then 

K-i  =  K;i^-a  K^^d  [  1  -H  fl  d + K;;^d  ]-i  d  (25) 


The  output  energy  of  adaptive  MFB  becomes 

1  +  a?d+  K-td, 


(Ea)  = 


'k'‘k 


(26) 


When  the  steering  vector  matches  the  signal  vector  (s  =  dt)  the  term  in  (..)  equals  zero,  giving 
the  following  equation  for  the  output  energy; 


{Ea)  =  + 


s+K-is 


(27) 


From  repeated  applications  of  Eq.  17  to  in  Eq.  21  (see  Appendix),  the  failowing  expression 
can  be  derived  for  the  inverse  of  the  total  noise  covariance  matrix: 


M(  a2d,d.+  1 

i_y 

j^k 


(28) 
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The  output  energy  of  the  adaptive  MFB  process  is  given  by 


(29) 


and  the  array  gain  is  given  by 


=  10  log 


*• 

• 

• 

■« 

N(\  +  aj) 

M 

1^2 
7  =  1 
jxk 

r  1 

7  7i 

1  +  JVa? 

L  7  J 

(30) 


Using  ^  (Eq.  42),  the  adaptive  MFB  output  energy,  assuming  all  interferences  are  at  the 
mean  sidelobe  level,  is  given  by 


{Ea) 


(31) 


and  the  array  gain  is  given  by 


AGj^  =  10  log 


N{\  +  aj) 


M 


i-Z 

7=1 


a? 


1  +^a? 

7 


(32) 


As  the  number  of  elements  increases,  N>M,  the  array  gain  increases  asymptotically  as 


AG^  "*  10  log  +  10  log  Oj 

N  >  M 


(33) 
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The  array  gain  is  then  bounded  by  10  log  N  plus  HNRj. 

The  array  gain  and  output  energy  from  adaptive  MFB  can  be  calculated  for  various  by 
using  Eq.  31  and  32,  and  by  assuming  a  distribution  of  interference  levels,  aj .  We  assume,  in 

flgure  2,  that  the  acoustic  field  consists  of  150  sources  of  varying  levels,  with  a  total  level  of 
79.95  dB.  In  this  example,  the  total  noise  level  is  80  dB.  The  total  input  interference  to  indepen¬ 
dent  noise,  IINRf,  is  assumed  to  be  10  dB. 


Figure  2.  Noise  field  assumption  for  figure  3.  Total  noise  level  =  80  dB. 
Independent  noise  level  =  69.95  dB.  150  interferences  uniformly  distributed 
in  energy,  totaling  79.95  dB. 


The  received  energy  when  the  array  is  focused  at  each  of  the  interferences  is  shown  in  fig¬ 
ure  3  for  various  number  of  array  elements,  N.  As  N  increases,  the  measured  noise  floor  drops 
and  more  signals  are  detected.  For  example,  with  1,000  hydrophones,  signals  approximately 
40  dB  below  the  80-dB  total  noise  level  can  be  detected,  achieving  a  40-dB  array  gain,  10  dB 
above  10  log  N. 

The  array  gain  in  this  noise  field,  calculated  by  using  Eq.  32,  is  shown  in  figure  4  versus 
N.  In  this  plot,  IINRt  is  varied  from  -30  to  +20  dB.  The  conditions  of  figure  2  and  figure  3, 
IINRj  ~  10  dB,  are  marked  with  “+”  in  the  figure.  From  figure  4,  measurements  with  1,000- 
element  volumetric  arrays  can  determine  the  IINR.  The  array  gain  achievable  with  a  larger  num¬ 
ber  of  elements  is  10  log  N  +  IINRt. 
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INTERFERENCE  (TARGET)  NUMBER 


Figure  3.  Received  level  from  adaptive  MFB  versus  number  of  elements,  N, 
when  the  processor  is  steered  at  the  interference. 


INTERFERENCE  (TARGET)  NUMBER 


Figure  4.  Array  gain  vs  number  of  elements,  N,  and  total-interference-to- 
independent-noise  ratio,  IIINRj. 
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PEAK  SIDELOBE  LEVEL  REQUIREMENTS  FOR  ADAPTIVE  MFB 

The  results  presented  above  are  expected  values  since  the  mean  sidelobe  level  of  -10  log  N 
was  used.  In  order  to  establish  array  design  criteria,  the  portion  of  the  noise  floor  energy  of  a 
given  interference  is  estimated  by  assuming  one  signal  (k),  one  interference  (i),  and  independent 
noise  in  Eq.  29.  This  is  given  by 


(£^>(st  d,-)  = 


1  - 


1 _ 

NajL^ 

1 


(34) 


The  second  term  is  the  contribution  to  the  noise  floor  from  the  interference.  It  is  plotted  in  fig¬ 
ure  5  versus  the  sidelobe  level  (Liy)  as  a  function  of  N  and  the  individual  interference-to-indepen- 
dent-noise  ratio,  IINRi  =  10  log  a? .  For  large  IINRi^  an  interference  contributes  1  dB  to  the 
noise  floor  for  peak  sidelobes  below  -6  dB  re  the  main  lobe.  A  few  sidelobes  as  high  as  -3  dB 
are  acceptable.  Arrays  can  be  designed  by  using  the  linear  matched-field  beampattem  statistics  in 
a  process  similar  to  random  arrays.^*^®4i,i2 


-20  -18  -16  -14  -12  -10  -8  -6  -4  -2 


SIDELOBE  LEVEL  (dB) 


Figure  5.  Contribution  to  the  measured  noise  floor  from  a  single  interference 
versus  the  Bartlett  sidelobe  level  at  the  interference  for  vaious^  X  IINRi. 


The  high  sidelobe  tolerance  of  adaptive  MFB  has  great  impact  on  the  design  of  arrays,  elimi¬ 
nating  the  need  for  large  numbers  of  accurately  positioned  high-tolerance  hydrophone  channels, 
which  are  required  to  obtain  the  low  sidelobe  levels  needed  to  detect  weak  signals  in  multiple 
interferences. 
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SENSITIVITY  OF  MATCHED-FIELD  BEAMFORMING  TO 

MISMATCH  ERRORS 

Errors  in  element  location,  propagation  modeling,  and  sound-speed  structure  measurement 
result  in  errors  in  the  prediction  of  the  signal  replica,  s.  These  mismatch  enors  cause  both  loss  of 
array  gain  and  bias  errors  in  estimating  the  source  position.^  The  latter  occur  when  a  near  p)erfect 
match  is  obtained  for  true  source  positions  that  differ  from  the  estimated  positions,  7,  ^  d.  The 
sensitivity  of  array  gain  and  output  spectrum  of  optimum  array  processors  has  been  thoroughly 
analyzed  by  Cox.^^  Though  derived  assuming  plane  wave  propagation,  his  results  are  directly 
applicable  to  matched-field  beamformers.  For  this  report,  we  examine  the  array  spectral  output 
degradation  (ASOD)  due  to  mismatch,  defined  as  the  ratio  of  the  energy  at  the  signal  peak  with 
mismatch  to  the  energy  with  no  mismatch: 


ASOD  =  10 


(35) 


This  definition  is  experimentally  estimated  using  a  strong  signal  at  a  known  position.  The 
energy  in  the  peak  nearest  the  known  location  is  measured  and  compared  to  the  predicted  energy 
with  no  mismatch.  The  definition  and  measurement  include  a  component  of  the  output  energy 
from  increased  sidelobes  caused  by  the  nulling  of  the  mismatched  signal.  This  section  describes 
the  array  gain  degradation  (AGD)  for  linear  and  adaptive  MFB  as  a  function  of  input  SNR.  It  is 
assumed  that  the  acoustic  field  is  generated  by  a  single-point  source  and  incoherent  noise.  The 
covariance  matrix  is  shown  in  Eq.  10. 

UNEAR  MFB  MISMATCH  EFFECTS 

In  practice,  the  signal  energy  measurement  includes  some  noise  energy,  as  shown  in  Eq.  12. 
Linear  anay  spectral  output  degradation,  A50D/,,  is  obtained  using  Eq.  11  and  12  in  Eq.  35: 


ASODi  =  10  log 


=  10  log 


oKs'^dd 

+s) 

+  J 

+ 

vC 

cos2(s,d) 

+ 

V' 

Na^ 

1  -1- 

X? 

Na^ 

(36) 


where  a*  = 


is  the  input  SNR  and 


cos^(s,  d) 


(s  +  d)(d  +  s) 
(s  +  s)(d+d) 


(37) 


is  the  generalized  cosine  between  the  steering  vector,  s,  and  the  signal  vector,  d.^^  The  latter  is  a 
measure  of  the  mismatch  signal  gain  degradation  which  results  from  nonalignment  between  the 
steering  vector  and  the  received  signal  vector.  It  is  also  the  ASOD  and  AGD  for  a  large  output 
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SNR.  The  linear  ASOD  for  various  output  SNRs,  Na\  is  plotted  in  figure  6  as  a  function  of 
mismatch  AGD,  10  log  cos^(s,  d).  We  have  assumed  equal  noise  on  all  of  the  elements,  A7’  =  i 
Note  that  the  output  SNR  must  be  20  dB  or  greater  for  a  direct  measurement  of  ASOD.  How¬ 
ever,  a  measurement  of  the  input  SNR,  X7*>  and  ASOD  can  be  used  with  Eq.  36  to  estimate  the 
mismatch  AGD  with  finite  SNR. 


Figure  6.  Bartlett  array  gain  degradation  versus  mismatch  array  gain 
degradation  for  various  output  SNRs. 


Following  the  work  of  Berman,^^  we  assume  that  s  and  d  differ  only  in  phase,  dj  =  sjer^K 
The  true  mismatch  is  given  by 


cos^(s,d)  = 


N 


7=1 


N 

1 

Lit=i 


N2 

N  N  N 

7*1  7=1*  *=1 


(38) 


N2 


If  it  is  assumed  that  the  mismatch  enors  are  represented  by  independent  Gaussian  distributed 
phase  errors  on  each  element  of  the  array,  given  by 


P(h  = 


1 


_1 

e  2*1 


(39) 
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where  is  the  RMS  phase  error  in  radians;  then 

=  e-^  (40) 

We  have  used  the  fact  that  the  probability  density  of  the  difference  of  two  independent  zero 
mean  normal  distributions  is  also  a  zero  mean  normal  distribution  whose  variance  is  the  sum  of 
the  component  variances. 

JV  N  N 

■*'  X  ^  ™  expected  value  of  the  output 

>-l  /  >-l  j-l  j*  k  k-i 

with  mismatch  is  given  by 


(cos^(s,d))  = 


1  +  (Nr-l)e-^ 


where  N,  =  is  the  effective  number  of  elements  illuminated  by  the  signal. 

y-i 

The  mismatch  AGD  from  phase  errors  is  plotted  in  figure  7.  It  can  be  seen  that  the  signal 
replica  at  a  range  near  the  true  range  must  be  calculated  to  ^  relative  accuracy.  This  implies  that 

relative  array  element  locations  be  known  to  better  than  ^  and  that  the  acoustic  models  and 
sound-speed  structure  measurements  generate  replicas  to  this  relative  accuracy. 


RMS  ERROR  {IN  WAVELENGTHS) 

Figure  7.  Expected  mismatch  array  gain  degradation  versus  RMS  phase 
error  in  wavelengths. 

Equation  41  can  also  be  used  to  derive  an  expression  for  the  mean  sidelobe  level.  Consider 
the  case  where  the  array  is  steered  away  from  a  source.  In  this  case,  the  mean  phase  error,  is 
large  and  the  mean  sidelobe  level  is  given  by 

(cos2(s,d))  =  r=  ^  (4: 
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ADAPTIVE  MFB  MISMATCH  EFFECTS 

Assuming  a  single  point  source  and  uniform  incoherent  noise,  the  covariance  matrix  is  given 
by  Eq.  10  with  A  =  I.  The  adaptive  MFB  array  spectral  output  degradation  due  to  mismatch — 
obtained  using  Eq.  9,  A-3,  and  38 — is  given  by 

ASODj^  =  10  log 

ASODj\  is  plotted  versus  mismatch  AGD  for  various  output  SNRs,  No?’,  in  figure  8.  At  high 
output  SNRs,  large  degradation  results  from  rather  small  mismatch.  In  this  case,  the  adaptive 
processor  is  nulling  the  portions  of  the  signal  that  do  not  match  the  replica.  However,  low-level 
signals  near  the  output  noise  level  (-10  log  N)  have  degradation  only  slightly  larger  than  the 
degradation  from  linear  MFB.  The  adaptive  processor  does  not  null  the  mismatch  portions  of 
low-level  signals. 


MISMATCH  AGD  (dB) 


Figure  8.  Adaptive  array  gain  degradation  versus  mismatch  anay  gain 
degradation  for  various  output  SNRs. 

Mismatch  errors  with  adaptive  MFB  thus  degrade  the  performance  more  than  linear  MFB 
only  at  high  SNR.  At  low  output  SNR,  the  mismatch  degradation  is  approximately  the  same  as 
with  linear  MFB.  The  dynamic  range  of  the  search  space  is  thus  compressed. 

The  degradation  at  high  SNR  can  be  reduced  by  bounding  the  magnitude  of  the  adaptive 
weight  vectors  using  the  white  noise  constraint  algorithm.^3>l^>^^>^^  These  methods  prevent  nul¬ 
ling  of  the  signal  in  the  main  beam  by  limiting  the  depth  of  the  nulls  in  the  focal  pattern.  For 
high-level  interferences,  the  contribution  to  the  noise  floor  shown  in  figure  5  will  be  greater  than 
shown,  depending  on  the  value  of  the  bound  of  the  weight  magnitude. 
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CONCLUSIONS 

Adaptive  mr^ched-field  beamforming  with  large  arrays  is  predicted  to  achieve  gain  equal  to 
the  directivity  index  plus  the  total  IINR.  An  example  of  a  physically  derived  distribution  of 
received  interference  levels  illustrates  that  measurements  with  volumetric  arrays  of  approxi¬ 
mately  1,000  hydrophones  can  determine  the  IINR,  enabling  an  estimation  of  achievable  array 
gain. 

Array  sidelobe  level  requirements  for  adaptive  MFB  are  derived,  showing  that  loud  interfer¬ 
ences  contribute  less  than  1  dB  to  the  noise  floor  for  peak  sidelobes  below  -6  dB  re  the  main 
lobe.  A  few  sidelobes  as  high  as  -3  dB  are  acceptable. 

Adaptive  matched-field  gain  degradation  from  replica  errors  is  shown  to  approach  the  linear 
matched-field  gain  degradation  at  low  output  SNR. 
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in  Eq.  21  is  of  the  form 
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c„  =  Xb^d;  +  I 

J  =  1 


(A-1) 


where  it  is  assumed  that  the  interferences  are  orthogonal: 


b/b*  = 


(A-2) 


This  appendix  shows  the  proof  by  induction  that 


w  b,b  + 
i_y  _i_i — 
1  +  6? 


(A-3) 


Define 


Co  +  bi  bj^ 
Cj  +  b2b^ 


(A-4) 


Cm  ~  +  bnibm 


Eq.  17  can  be  used  to  show  that 


*  l+b;c-‘,b. 


(A-5) 
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Thus,  using  Eq.  A-2 
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=  I  + 


1  +  b\ 
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(A-6) 


If  we  assume  that 


then,  using  Eq.  A-5 


m-\ 


b,b.+ 

I-  y  — — — 
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Y  b^b  + 


m-1 
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U'l  J 
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All  the  sums  in  the  third  term  are  zero  by  the  orthogonality  assumption,  Eq.  A-2,  yielding 
Eq.A-3. 
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